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We present the most recent lattice results for the lowest-order hadronic contribution to the muon 
anomalous magnetic moment using 2+1 flavor improved staggered fermions. A precise fit to 
the \ow-q 2 region of the vacuum polarization is necessary to accurately extract the muon g — 2. 
To obtain this fit, we use staggered chiral perturbation theory with the inclusion of the vector 
particles as resonances, to evaluate the vacuum polarization. We discuss the preliminary fit results 
and attendant systematic uncertainties, paying particular attention to the relative contributions of 
the pions and vector mesons. 
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Figure 1: The lowest order hadronic contribution to the muon g—2. The blob represents all QCD diagrams. 



Currently, experimental measurements of the muon g—2 have reached a level of high precision 
[jl|]. As a check on the current standard model of particle physics, a competitively precise theoretical 
estimate for this quantity is needed. The current theoretical calculations come from dispersion 
relations relating the photon vacuum polarization to the cross-section for either e + e~ -^hadrons or, 
by using isospin symmetry, T decay to hadrons. The agreement between the theoretical calculations 
coming from e + e~ data to experiment is within 2.7 standard deviations, and this is lowered to 
1.4 standard deviations if one includes the results incorporating z decay [Q, ||]. The need for 
a determination completely from first principles (i.e., using the lattice) is desirable to check the 
current discrepancy between theory and experiment. 

We can extract the muon anomalous magnetic moment by calculating the full muon-photon 
vertex, where the correction to the tree-level result is given by = (g — 2)/2 = F2(q 2 = 0), where 
F2(q 2 ) is the form factor proportional to [y^,/ v ] in the full vertex, and q is the momentum trans- 
ferred to the photon. This can be calculated in perturbation theory in the electroweak theory, with 
the first contribution coming in at 0(a), where a = e 2 /An is the fine-structure constant. 

The first effects from QCD come in at 0(a 2 ) and are shown in Fig. |l|. The hadronic contri- 
butions are 7 x 10~ 5 smaller than the leading corrections and account for most of the discrepancy 
with experiment. Due to the strength of the sttong coupling constant, these terms cannot be reliably 
calculated in perturbation theory at the relevant energies, so other methods must be employed to 
calculate them, and this is precisely where the lattice can be useful. Currently, the next hadronic 
contribution, coming from light-by-light scattering at 0(a 3 ), is also being studied using lattice 
techniques [ffl|. 

Lattice QCD provides a first principles method to calculate the hadronic contribution to the 
muon g—2. We can calculate the QCD contribution to the muon g—2 coming from the vacuum 
polarization of the photon using Refs. \K |6l [7||, 

a^ d =(^) 2 J~dK 2 f(K 2 )n(K 2 ), (1) 

where Yl(K 2 ) is the vacuum polarization of the photon, which we calculate on the lattice, and f(K 2 ) 
is a kernel given in Ref. [Q]. Since f(K 2 ) diverges as K 2 — > 0, we must be able to calculate the low 
momentum region of the vacuum polarization very precisely, and thus large lattices will be needed 
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Figure 2: The two diagrams that contribute to the photon vacuum polarization at one loop in SZPT. 

here. This is also why perturbation theory is not reliable, since the largest contribution to aj^ had 
comes from the low-energy regime. 

A detailed discussion of the lattice calculation can be found in Refs. (EL EJ. We used lattices 
generated by the MILC collaboration using 2+1 flavors of light dynamical improved staggered 
quarks. The results shown here are those coming from three of the "fine" lattices (a 0.086 fm) 
with a lattice strange sea quark mass of 0.031. The light sea quark mass in lattice units is 0.0062 
on the lattice volume of 28 3 x96 and 0.0031 on the 40 3 x 96 volume. 

In order to fit to the low-momentum region, we use Staggered Chiral Perturbation Theory 
(SZPT) [®J coupled to photons. The generalization of the standard staggered chiral Lagrangian 
to include external photons amounts to changing the derivative operator to a covariant derivative: 
d^L —>■ D^L = d^L + ieA^ [Q, E] , where £ is the 12 x 12 matrix of pseudo-Goldstone bosons, Q is 
the light quark charge matrix in units of e and is the photon field. Recall the chiral symmetry of 
SZPT is an SU (12) l xSU (12) r symmetry for 3 flavors of staggered quarks, and under this symme- 
try £ transforms as £ — > LLR^ with L G SU (12) l and R G SU (12) r. Our (Euclidean) Lagrangian in 
this case is given by 

f- uf 2 
JSf = y Tr^ZD^r 1 "] - ^-Tr [JTL + L^}+a 2 Y (2) 

where Y is the taste-symmetry breaking potential for multiple flavors of staggered quarks [B] and 
j$ is the light quark mass matrix. No additional terms arise at this order in the Lagrangian which 
violate the taste symmetry while being coupled to the external photons. Additionally, we drop the 
additional singlet mass term from Ref. [||] as it will not be relevant for the current calculation. 
Using this, we can calculate the one-loop contribution (Fig. which is given by 

where x = 4m 2 / p 2 , i £ {n + ,K + }, and the sum over t is over the 16 tastes of mesons. The masses 
for these mesons in terms of the Lagrangian parameters are given in Ref. [Jsj] and their measured 
values on the lattice can be found in Ref. This is the only place where the effects of taste 
violations are seen. There are no neutral mesons in the loops, so we do not get effects from the 
two-point hairpin insertions [Bj. Additionally, aside from an overall constant, A, there are no free 
parameters in this expression, as all the meson masses have been determined in lattice simulations 

Hi- 
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Figure 3: The diagrams that contribute to the photon vacuum polarization with the p in SZPT at (a) tree 
level and (b) one loop. The double line is the vector mesons and the single line is a pion or kaon as before. 



The problem, however, is that this is two orders of magnitude too small compared with the lat- 
tice data. This indicates that there must be other effects which dominate the low energy physics of 
the hadronic contribution to the vacuum polarization. We take a lesson from quenched simulations 
10, where the vacuum polarization can be completely described by the p. In that case, there is a 
term coming from the vector bound state and an additional term coming from the two-particle con- 
tinuum. This implies that perhaps in the dynamical case, there may be large effects coming from 
vector mesons, despite their large masses, which are not included in standard chiral perturbation 
theory. 

The model we use here is that of the resonance formalism of Ref. JTT|], where the relevant 
interaction Lagrangian is 

ifvec = ^Tr [y M v(aF MV a t + a^a)] . (4) 

fv is the tree-level vector decay constant, Vn V is the 12 x 12 matrix of the 9 vector mesons, a 2 = £, 
and F^ v = eQ(d ll A v — dyA^). Although we denote the vector fields using an antisymmetric tensor 
field, it can be shown that there are still the correct number of degrees of freedom for each particle 

EDl- 
Under the chiral SU(12)l x SU (12) r symmetry, we have 

Vftv -» , o -» LoU 1 = UoR 1 (5) 

where U is a spacetime-dependent SU (12) matrix defined by this equation. 

The leading order contribution here is the tree-level diagram shown in Fig. |||(a), and is given 

by 



n v (p 2 



a (4nf\ 



v, 



\2 



An 



3 1 

+ ■ 



p 2 + m 2 p0 p z + m z a 



(6) 



We have not included the taste index on the vectors, since empirically the taste-breaking for the 
vectors is negligible [g]. Again, we see that if we measured fy and the vector masses on the lattice, 
there are no free parameters in this expression. Also, although the masses of the vectors are much 
heavier than the pion, there is the enhancement factor of (4nfv) 2 in the numerator here, which we 
can see makes this the dominating factor in the overall result. Ref. [12] estimates fy to be about 
200 MeV, so this term is roughly 0(1) in the chiral power counting scheme, and thus will dominate 
over the one-loop pion term. 
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Figure 4: Fits to the lattice data using the tree-level vector formula for the three different quark masses for 
the \ow-q 2 region. 



Since we worked to one-loop order in the pion sector, we should do so in the p sector, as 
these should be roughly the same order. The only contributions are tadpole corrections to the p — y 
vertex, Fig. ^(b). We get 



p 2 + m 2 , p 2 + rrA p 2 + in 



(7) 



We can see again that the taste violations here enter quite trivially as they did in the pure pion 
sector. 

Figure || shows the data for the three different masses used with the results from fitting to 
Eq. (g). Although these fits only include the tree-level p results for now, the inclusion of the one- 
loop formulae changes the fit only slightly. We can see that Eq. (^) fits quite well, except for lighter 
quark masses and the smallest value of q 2 , where the fit function undershoots the data. There are 
several reasons for why this may occur. We do not include the Naik term when calculating the quark 
propagators, though it is hard to see how this would affect the low q 2 , long-distance, behavior. It is 
possible the effect is due to finite volume, though in the free theory the vacuum polarization turns 
over the other way from small q 2 . Finally, it may simply be due to poor statistics. 

Using the results of these fits, we find when extrapolating to the physical values of the quark 

9 had VP 

masses that the hadronic contribution to the muon g—2 (from the q < 1 GeV regime) is ' = 
657(20) x 10~ 1() . The perturbative contribution coming from the q 2 > 1 GeV regime is a ^ ad V P,pert < 
10 x 10~ 10 , and the final result is the sum of these two terms. The error quoted here is statis- 
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tical only. This is somewhat lower than the results from dispersion relations, which is ' = 
693.4(5.3) (3.5) x 1(T 10 . The difference with the previous lattice result of a^ yP = 545(65) x 1CT 10 
[j7|] comes mainly from the inclusion of the 0.0031 light quark mass result which allows an extrap- 
olation, using Eqs. (|J) and to the physical quark mass. This accentuates the fact that a theo- 
retically motivated fitting function is necessary to obtain an accurate determination of the low-g 2 
region. We also note that the final value of the lowest order hadronic correction may increase fur- 
ther if the observed overshoot of the fit function by the lattice data turns out to be physical. Finally, 
we caution that the disconnected diagrams (connected by gluons) that contribute to the lowest order 
vacuum polarization have not been included yet. These are Zweig suppressed in general and vanish 
exactly in the SU (3) flavor limit, but may contribute on the one-to-few percent level. 

This work was supported by the U.S. DOE. We thank NERSC for the computational resources 
expended on this project. TB thanks Michael Ramsey Musolf for helpful discussions. 
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